Determination of intrinsic switching field distributions in perpendicular recording 
media: numerical study of the AH(M, AM) method 
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We present a numerical study of the AH{M, AM) method and its ability to accurately determine 
intrinsic switching field distributions in interacting granular magnetic materials such as perpendic- 
ular recording media. In particular, we study how this methodology fails for large ferromagnetic 
inter-granular interactions, at which point the associated strongly correlated magnetization rever- 
sal cannot be properly represented by the mean-field approximation, upon which the AH{M, AM) 
method is based. In this study, we use a 2-dimensional array of symmetric hysterons that have 
an intrinsic switching field distribution of standard deviation a and ferromagnetic nearest-neighbor 
interactions J. We find the AH{M, AM) method to be very accurate for small J/a values, while 
substantial errors develop once the effective exchange field becomes comparable with a, corroborat- 
ing earlier results from micromagnetic simulations. We furthermore demonstrate that this failure 
is correlated with deviations from data set redundancy, which is a key property of the mean-field 
approximation. Thus, the AH{M, AM) method fails in a well defined and quantifiable manner that 
can be easily assessed from the data sets alone. 

PACS numbers: 75.50.Ss, 75.60.Ej, 75.75.-(-a 
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I. INTRODUCTION 



One of the key challenges in advancing the nanotech- 
nology of magnetic recording is the optimization of 
recording media and its physical properties^. This chal- 
lenge is particularly demanding because magnetic record- 
ing is by its nature a local process. Thus, it is not so much 
the average physical properties that are crucial, but the 
distributions of such properties that determine contin- 
ued technological advancement and success^. In general, 
it is important to devise recording media structures that 
have very homogeneous properties on recording relevant 
length scales, so that a position independent physical de- 
scription of all properties and magnetization processes 
is appropriate. However, this can only be achieved to 
a limited degree and it is therefore essential to have ex- 
act knowledge of the corresponding parameter distribu- 
tions. One of the most crucial properties is the intrinsic 
switching- field distribution D{Hs) of the media grains 
because it defines the recording quality of a media layer 
in both magnetic stability and the achievable recording 
density^. Hereby, one has to realize that it is not the 
macroscopic switching field distribution D^{Hs) in a uni- 
formly applied field that is relevant, but the local distri- 
bution D{Hs) of switching fields in a recording process, 
which takes place in a narrowly defined field geometry. 
The difference between D^{Hs) and D{Hs) is caused by 
the inter-granular interactions between the media grains. 
In particular, for perpendicular media the dipolar inter- 
action is large and dominates the behavior of hysteresis 
loops M{H). Thus, the knowledge of D„^{Hs), which 
can be derived from the slope of M{H), is insufficient for 
recording performance predictions. 



Over the years, several methodologies have been devel- 
oped to determine D{Hs) with varying success^i^i^i^i^iiSii^. 
Most of these methods apply a measurement scheme, in 
which magnetization reversal of media grains is probed 
starting from different magnetization states to allow for a 
variation of the grain-to-grain interaction under measure- 
ment conditions. Such procedures should then in princi- 
ple and under certain limiting conditions allow for a sep- 
aration of the intrinsic switching field distribution and 
the inter-granular interactions. 

In this paper, we study the reliability of the recently 
developed AH{M, AM) method, which has been used 
successfully in analyzing progress in recording media fab- 
ricationli^. The method itself is a generalization of an 
earlier measurement technique, the AHq method?, but 
has the advantage that it allows the determination of the 
entire D{Hs) distribution and not just a single character- 
istic parameter. Furthermore, it enables oversampling, 
which makes consistency checks feasible and gives one 
the opportunity to quantify the confidence level of mea- 
surement results. In contrast to the also quite robust 
method developed by van de Veerdonk et al. — i^, it does 
not rely on a specific distribution form or the rather lim- 
iting assumption that interactions can be removed from 
the problem by a simple de-shearing of the major hystere- 
sis loop. Recent data indicate that this very assumption 
appears to be an overly simplistic view of inter-granular 
interactions^. Another method that has recently gained 
certain popularity is the FORC-method'*'^, which is very 
similar to the measurement of Preisach distributions^-. 
However, this technique cannot really be compared to 
the previously mentioned methods, because it represents 
for the most part a data transformation tool and does 



2 



not appear to allow a self-consistent way of extracting 
microscopic information such as D(Hs), because all inter- 
actions are removed from the model in an ad hoc fashion 
simply by definitional^. 

Here, we present a numerical study of the 
AH(M, AM) method at zero temperature. We 
show that even though this method approximates the 
inter-granular interactions on the mean-field level, it can 
predict its own reliability correctly. The paper is orga- 
nized as follows. In Sec. [TTl we give a brief introduction 
to the AH{M, AM) method. In Sec. |TTT1 we discuss the 
model used to numerically calculate the magnetization 
curves M{H). In Sec. IIVI we introduce reliability 
measures to quantify the failure of the AH{M, AM) 
method. In Sec.|Vl we present and analyze our numerical 
results and a summary is given in Sec. IVII 



II. 



AH{M,AM) METHOD 



The AH{M,AM) method assumes that the effec- 
tive field at each grain can be written as H^s = 
H{M) + Hint{M) with H{M) being the external field and 
Hint{M) the volume-averaged interaction field at magne- 
tization M . Assuming that M is normalized to its satu- 
ration value Afg, the upper branch of the major hysteresis 
loop driven by the major loop external field H, and the 
recoil curve starting at Micv = 1 — AM driven by the 
recoil curve external field Hi are given as: 



/.-[//(Af)+ffi„t(M)] 

M = 1 - 2 / D{Hs) dHs (la) 

^ — OC 

M =1- AM - 2 / D{Hs) dHs (lb) 

respectively, i.e. field integrals over the intrinsic switch- 
ing field distribution D{Hs). Defining an integral func- 
tion I{x) = D{Hs) dHs, one finds 



H{M) 



1 - M 



Hint{M) 



(2a) 



Hi.iM)-^I-'{^^^^^)^HUM) (2b) 

with being the inverse of the integral function. 
Therefore, one can derive 



AH{M,AM) = H,{M)-H{M) 
1 - M 



r 



l-M -AM\ 



(3) 



as a closed functional form for the field axis distance AH 
between major and recoil curves, a quantity that is illus- 
trated in Fig. [3) Important is hereby that within the 
framework of this approach AH is independent from the 
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FIG. 1: Gaussian, Lorentzian and Lognormal D{Hs) distribu- 
tions with the same disorder parameter a. For both Gaussian 
and Lognormal distribution, the standard deviation is set to 
be 4.0 and the mean value 20.0. For the Lorentzian distri- 
bution, we set the width at half-maximum to be 4.0 and the 
center of the distribution to be 20.0. 



grain interaction, which allows for a direct experimental 
access to D{Hs)- For certain parameterized distribution 
functions, one can derive analytic expressions for AH. 
For example, for the Gaussian, Lorentzian and Lognor- 
mal D{Hs) distributions (see Fig. [1]), 
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we find 



AHg{M, AM) = V2a (erf"^(M + AM) - erf"^(M)) 

(5a) 

AHl{M, AM) = ^ |tan[|(A/ + AM)] - tan[|A/]| 

(5b) 

Ai7LN(Af, AM) = exp [fl - V2deTr\M)] 

- exp [/i - V2deTr\M + AM)] 

(5c) 

Results for the Gaussian and Lorentzian distributions 
were reported previously^. Here, we introduce the dis- 
order parameter a for a general D(Hs) distribution, a is 
defined to be the standard deviation of a given distribu- 
tion, such as Gaussian and Lognormal. However, for the 
Lorentzian distribution neither variance nor higher mo- 
ments are defined, so that we need to quantify the disor- 
der level in another form. For the Lorentzian, we define 
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the disorder parameter to be the distribution width at 
the half- maximum and the mean to be the center of the 
distribution. For the distributions given by Eq. [4al [4bl 
and del we then have the disorder parameter: 
cr^ = w, a^^ = e'^+*'/2(e*' - 1)^/^ and the mean value 
= ho, h}^ = ho, h};'' = e'^+*V2. 

By making a least-squares fit of the AH{M, AM) 
curves to the above formulas, one can extract the key 
features of D{Hs). Note that both AHg{M,AM) and 
AHi^{M, AM) have no ho dependence. But for the Log- 
normal distribution, AHi,^{M, AM) depends on both jl 
and a, and therefore it has both a and ho dependencies. 

It is easy to prove that simply shifting a general distri- 
bution will not alter the AH data. For this, we consider 
the case of a general distribution D{Hs): if we shift it 
towards right by an amount Ho, then the integral func- 
tion I{x) of the new distribution D{Hs) — D{Hs — Ho) 
is given by I{x) — I{x — Ho) and its inverse I~^{y) — 
Ho+I~-^{y). It is then clear that this shift will not change 
the AH{M, AM) formula at all, since the Ho terms will 
cancel according to Eq. [3] 

Also, we find that both AHg{M,AM) and 
AHi^{M,AM) are symmetric with respect to 
M = -AM/2 while AHi^^{M,AM) is not sym- 
metric. This is consistent with the original distri- 
bution: Both Dg{Hs) and Di,{Hs) are symmetric 
around Hq while I?ln(^^s) is not. Actually, for 
any D(Hs) distribution being symmetric about ho, 
one finds that /(ft-o + x) = 1 — /(/lo ~ x) and 
/~^(l/2 + y) = -/-i(l/2 - y), from which it is 
easy to prove that AH[M, AM) is symmetric about 
M = -AM/2, i.e. 



AH{ — + M', AM) = AH{ M', AM) (6) 



III. HYSTERON MODEL 

For our numerical studies, we model each grain as a 
hysteron, which is the simplest mathematical construc- 
tion for the description of a hysteretic system. Each hys- 
teron will generate a rectangular hysteresis loop in an 
applied field H as shown in Fig. [21 We assume the tran- 
sition of each hysteron is infinitely sharp and it has no 
additional field dependence, such as finite susceptibility, 
for instance. The half-width of the rectangular hystere- 
sis loop is referred to as the intrinsic switching field of 
the hysteron, which is a well defined property of each 
individual hysteron. 

We also assume that the hysteron switching fields 
are symmetric around zero field, i.e. there is no bias. 
A hysteron with zero bias is called a symmetric hys- 
teron and is consistent with the time reversal symme- 
try that ferromagnets exhibit^i. These hysteron prop- 
erties can be regarded as fairly good representation of 
perpendicular media grains, because they have relatively 
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FIG. 2: A symmetric hysteron S with intrinsic switching field 

high magnetic anisotropy and the applied field in typical 
characterization measurements is applied along the easy 



We further assume that there is no time structure to 
the hysteron switch itself. This should be an appropriate 
picture as long as one considers field change rates that 
are much slower than single grain reversal times. Given 
that such reversal times are typically of the order of sev- 
eral hundred pico seconds-^, this condition is generally 
fulfilled in typical measurement setups^'^. 

In our model the ferromagnetic layer system is then 
represented by a simple square lattice of symmetric hys- 
terons with periodic boundary conditions. Note that a 
square lattice is not necessarily a very good approxima- 
tion of actual media structures, in which grains typically 
have coordination numbers of 5 or 6^. However, we ne- 
glect this detail of actual media structures for reasons of 
simplicity and furthermore assume that deviations from 
a more realistic grain structure are not fundamentally 
altering the overall significance of our study. Under the 
assumptions that hysterons {Si = ±1) have an intrinsic 
switching field distribution D{Hsi) with > 0; inter- 
act ferromagnetically with their nearest neighbors with 
strength J and experience a uniform external field H , the 
Hamiltonian of the system can be written as 

n = -Y, JS^SJ-Y,(H + sgn{S,)Hs,) (7) 

<i,j> i 

We denote this as the random switching field hysteron 
model (RSFHM). Note that the RSFHM is very similar 
to the random field Ising model (RFIM): 

n^-J2 J S'^sj -J2iH + lM) (8) 

<i,j> i 

where the spins Si = ±1 interact ferromagnetically with 
their nearest neighbors with strength J and experience 
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a uniform external field H and a local quenched field 
hi- A simple mapping: sgn(5'i)iJsi enables us to 

calculate M{H) curves of symmetric hysterons with the 
developed algorithms used for previous RFIM computa- 
tional worki^ii^. Details of the close relation between 
these models are described in Appendix [XI 

In this paper, we evaluate different types of D{Hs) 
to study the reUability of the AH{M, AM) method. In 
our simulation, we set the ferromagnetic nearest-neighbor 
coupling strength J = 1 and tune the disorder parame- 
ter a. Since J = 1, tuning a is equivalent to tuning 
a/ J. In this sense, a small (big) a corresponds to strong 
(weak) nearest-neighbor interactions of hysterons. We 
calculate the M{H) curves, both major hysteresis loop 
and recoil curves, for system sizes up to 1000^ and a 
values from 1.6 to 1000. Note that we are ignoring dipo- 
lar interactions in our model (Eq. [7]), even though we 
know that they are substantial in real structures. But 
a previous micromagnetic study demonstrated that the 
dipolar interactions can very well be treated within the 
mean-field AH{M, AM) method and do not cause any 
significant precision problems, while exchange interac- 
tions did^. Thus, we focus here on the effect of inter- 
granular exchange interactions only, because they repre- 
sent the much more serious problem for the reliability of 
the Ai/(M, AM) method. Also, we are acutely aware of 
the fact that the assumption of a uniform exchange cou- 
pling constant J in Eq. [7] is a substantial simplification 
of the problem, if one compares our model to real mate- 
rials structures. However, a lateral varying J would be 
an extension of the present model and will be the topic 
of future work. 



IV. RELIABILITY MEASURES 

In numerical test simulations of characterization meth- 
ods, one can introduce a simple type of reliability mea- 
sure, namely the deviation of the retrieved fit parameters 
from the exact solution given by the input parameters of 
the simulation. However, for experimental data such a 
reliability measure is generally not available simply be- 
cause the materials parameter are not known, but are 
actually what one wants to extract by means of this very 
experimental method. Therefore, in a typical experimen- 
tal situation one can only guess the overall measurement 
reliability by identifying parameter ranges, in which ap- 
proximate methods are expected to work with a certain 
level of precision, which may or may not be clearly iden- 
tifiable. One broadly used method to identify the relia- 
bility of extracted materials parameter is the similarity, 
with which a fit that is based on a certain physical model 
can reproduce the experimental data. 

As we will demonstrate in the following, we en- 
counter an unusually good situation for the AH{M, AM) 
method, in that a consistency test of the data set alone 
can be used to validate the applicability of the method 
and estimate its overall reliability. We will furthermore 



show that this consistency test does not depend on any 
specific distribution shape but just the mean-field ap- 
proximation in general, i.e. the underlying assumption 
of the AH{M, AM) method. 

In this section, we discuss the reliability measures of 
the AH(M, AM) method for an arbitrary type of D{Ha,). 
The reliability is characterized by two types of measures: 
(1) conventional quality measures for numerical fits such 
as the percentage difference between the fitting and the 
actual parameter (Pd) and the square of the multiple 
correlation coefficient (i?^); (2) the average deviation 
from redundancy measure (r), which represents the above 
mentioned self-consistency test. 

A. Fit Quality 

In simulations, one knows the input parameters, in our 
case the input D(Hs) distribution, exactly. Furthermore, 
we have also derived AH{M, AM) formulas for certain 
specific D{Hs) distributions as shown in Sec. |TT1 Thus, 
by fitting these analytical formulas to the AH{M, AM) 
curves obtained from our 2D-RSFHM hysteresis loop 
simulations, we can get an estimate of the input D(Hs) 
distribution. Obviously, if the input D{Hs) is recovered 
by the fitting procedure with high accuracy, then the 
AH{M, AM) method works. To quantify the reliability 
of the AH{M, AM) method, we introduce the following 
fit quality measures. The most important fit quality mea- 
sure, denoted as Pd, is the percentage difference between 
the parameters obtained from a least-square fit and the 
actual input parameters into our simulation. It describes 
how well the AH{M, AM) method can indeed retrieve 
the information sought. Focusing on the disorder pa- 
rameter of D{Hs) as the most crucial fit parameter, we 
define 

Pd - [at - a)/a (9) 

Here CTf is the fit value of the D{Hs) distribution disorder 
parameter as defined in Sec. [Ill while a is the input value 
of the same parameter. 

Another fit quality measure that can be utilized here 
is the square of the multiple correlation coefficient R^, 
which measures how successful a fit and fit function is in 
explaining the datai^. It is defined as 



Yi and Yi are hereby the simulation result and the fit 
function value for AH{M, AM) at a given data point 
(Af,;,AMi), respectively. Y is the average value of Yi. 
According to the definition of P^, we know that as B? 
approaches 1, the fit is a better and better representation 
of the data set. Hereby, it is important to emphasize two 
points in calculating Pd and E? . First, we are comparing 
a mean-field theory to numerical simulations that contain 
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the complexity of magnetization reversal in its complete 
detail. Second, we are using a finite size grid in our sim- 
ulation while the analytic theory is derived for infinite 
systems. Naturally, the finite size will affect the E? and 
P<j calculation. Particularly, for data points near the be- 
ginning or end of the reversal curve, where only a few 
grains (hysterons) are reversed, the analytical theory for 
infinite hysteron numbers might not be accurate at all 
for the description of a finite system, independent from 
the validity of the mean-field approximation. 



B. Deviation From Redundancy 

As already mentioned above, fit quality measures such 
as E? are not necessarily foolproof because physical mod- 
els and corresponding fit functions may be used in a 
regime, for which the underlying theory does not apply 
anymore. In such cases, data fits and extracted materials 
parameters might appear very accurate, while they are 
not. Thus, it would be a tremendous help, if an indepen- 
dent data set evaluation would be available that allows 
a separate measure of the suitability of the underlying 
theory. Specifically here, this evaluation should tell us 
how good an approximation the mean-field theory is for 
any given data set, so that we know to which confidence 
level we can rely on the AH{M, AM) method. 
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We find that it is indeed possible to define such a quan- 
tity. To do so one has to recognize that within the mean- 
field approximation AH{M, AM) data sets contain re- 
dundancy. The data set redundancy can be seen from 
Fig. [3l For illustration purposes, we pick six points: 
(A and P), (B and U), (Q and V) located on the ma- 
jor loop, the i-th recoil curve and the j-th recoil curve, 
respectively. The i-th and j-th recoil curve start at 
-^^icv = 1 — AAIi and 1 — AMj , respectively. Without lim- 
iting the generality of this consideration, we furthermore 
assume that AMi < AMj and require that the six points 
satisfy: Ma = Mb = M, Mp = Mq= M- AMj + AMi 
and Mtj — My = M — AMj where M is an arbitrary 
value within f— 1, 1 — AMA. We then define 



AH,iM)^HB-HA (11a) 
AHj{M - AMj + AMi) ^Hq - Hp (lib) 
AHj{M - AMj) - AHi{M - AM,) =Hv - Hu (11c) 

Within the mean-field approximation, it is easy to 
prove that 

[Hb - Ha) + {Hy - Hij) - (iJq -Hp) = (12) 
as shown in Appendix [B) More generally, one finds that 



AH,{M) + AHj{M - AMj) - AH,{M - AMj) - AHj{M - AM^ + AM,) = 



(13) 



This data set redundancy is due to the fact that successive recoil curves are not fully independent and contain 
repeated information. However, Eq. [T3] is derived under the assumption of the mean-field approximation and is only 
true if the mean-field approximation is indeed fulfilled by the data set. For general data sets this property is not 
conserved. Therefore, we can define an M-dependent measure as 



,{M)^ 



AH,{M) + AHj(M - AMj) - AH,{M - AMj) - AHj{M - AM^ + AM,) 
AH,{M) + AHj{M - AMj) + AH,{M - AM,) + AHj{M - AMj + AM,) 



(14) 



r 



that monitors deviations from the mean-field approxima- 
tion based upon the above redundancy criterion (Eq. llSp . 
Eq.[Tlhas an M definition range of [AM^ - 1, 1 - AMi]. 
For a general set of multiple recoil curves, the average 
deviation from redundancy measure can be defined as 



Note that in calculating r, we are comparing the data 
set with itself. Thus, the validity of the underlying mean- 
field approximation used for all the AH{M, AM) data 
fits, can be assessed independently and from the data set 
alone. 



(15) 



with n being the total number of all the possible {i,j) 
pairsii. Thus, r is a quantitative measure that allows an 
accurate and independent check of how close or far any 
AH{M, AM) data set is from fulfilling the mean-field 
approximation. 



An additional benefit in the calculation of r is that the 
finite size inaccuracies at the definition range boundaries 
for the recoil curves will cancel out, at least to some 
degree. See Eq. [H For M = 1 - AMi, the finite size 



inaccuracies in AHi (M) and AHj (M - 
cancel out. Similarly, for M = AMj 



- AMj + AMi) will 

- 1, the finite size 



inaccuracies in AHj{M-AMj) and AHi{M-AMj) will 
cancel out. In this sense, the deviation from redundancy 
measure is more robust than the fit quality measures. 
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up to 1000^. For a Lorentzian, we choose /iq/c = 5 x 10^ 
for systems of size 100^ and /iq/c ~ 2 x 10^ for sys- 
tems of size 1000^. To avoid the long negative tails of 
the Lorentzian, we can truncate the distribution instead 
of choosing a huge ho /a ratio, creating a new type of 
distribution, which we refer to as truncated Lorentzian 
distribution Di^^{Hs) in the following. Further details, 
including the definition of a, for this distribution are de- 
scribed in Appendix [Cl For this Di^^{Hs) distribution, 
we also choose ho /a = 5 in our calculations. For the 
Lognormal distribution, there is by definition no distri- 
bution density for negative fields. But we still choose 
/lo/f = 5 to make it comparable with the Gaussian and 
the truncated Lorentzian distribution. 



A. Comparison with the mean-field approximation 



AH 




FIG. 3: (Top) The major hysteresis loop and three recoil 
curves. Throughout the paper, M is normalized to its sat- 
uration value Ms and H is normalized to the coercive field 
He- The first two recoil curves start at Mrcv = 1 — AMi 
and 1 - AMj , respectively, with AMi < AMj . For the six 
points (A and P), (B and U), (Q and V) shown in the figure 
with Ma = Mb = M, Mp = Mq = M - AMj + AMi and 
Mu = My — M — AMj where M is an arbitrary value within 
[—1, 1 — AMi], one can prove the existence of data redundancy, 
i.e. the equality {Hb ~ Ha) + {Hv - Hij) = {Hq - Hp) 
within the mean-field approximation. (Bottom) The corre- 
sponding AH {AI, AM) data for the three recoil curves. 



V. SIMULATION RESULTS 

To show that the AH{M,AM) method fails re- 
producibly in a well-defined manner, we calculate 
AH{M, AM) data sets for 4 different distributions. 

For all these distributions, we tune the disorder pa- 
rameter a from 1.6 to 1000, but keep the ratio ho /a to 
be a positive constant. This is done to avoid any nega- 
tive switching fields, which would otherwise describe non- 
physical behavior in violation of the second law of ther- 
modynamics^-. For a Gaussian, hg/a = 5 is generally 
big enough to avoid any negative Hg for system sizes of 



Key results of our numerical hysteresis loop calcula- 
tions for all these different switching field distributions 
are shown in Figs. [3]- [71 In each case, we show plots for 
a = 1.6,5 and 50 only to illustrate the general trends. 
All the calculations shown here are done in _D = 2 di- 
mensions for linear system size L = 1000, i.e. — 10^ 
hysterons. 

At first, we discuss the results for the Gaussian D{Hs) 
distribution in detail. Fig. [3] displays the results for dif- 
ferent (t's in different rows: (Top) cr ~ 1.6. (Middle) 
a — 5. (Bottom) cr = 50. For each a, we calculate a 
complete set of M{H) curves, both the saturation hys- 
teresis loop and recoil curves, as shown in the left column 
of Fig. [31 Note that here and throughout the paper, M 
(or AM) is normalized to the saturation value Ms and 
H (or AH) is normalized to the coercive field Hq- In 
particular, we choose 5 equally-spaced recoil curves, for 
which the distance to saturation is given by AM^ — i/3. 
From the left column of Fig. [H one can see that the hys- 
teresis curves get broader for higher cr. Note that it is 
not the larger cr itself that causes this effect, because this 
type of broadening is taken out due to the normaliza- 
tion of H with He, and the constant ho /a ratioi^. The 
difference in shape here actually reflects the fact that 
for lower cr, one gets correlated magnetization reversal 
which sharpens the macroscopic switching field distribu- 
tion substantially. 

In the middle column of Fig. [51 we show the corre- 
sponding AH{M, AM) curves (solid lines) derived from 
the simulated M{H) curves, as well as the mean-field 
approximation of the AH{M, AM) curves (dotted lines) 
calculated from Eq. I5al The mean-field curves are cal- 
culated by using the exact input parameter and are not 
a least square fit. This allows for a clear illustration of 
the deviations from mean-field behavior. From Fig.[3I^d), 
we see that for small a (corresponding to strong hys- 
teron interactions) the difference between the numerical 
result and the mean-field approximation is large. For 
intermediate a (corresponding to intermediate hysteron 
interactions), the difference diminishes but is still visi- 
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ble, especially near the negative saturation M = — 1, as 
shown in Fig. [D^e). For high a (corresponding to weak 
hysteron interactions), the difference is so small that it is 
not visible in Fig.lD^f). It should be mentioned that due 
to the constant /iq/c ratio and the normalization of AiJ, 
the mean-field Ai?(M, AM) curves look almost identi- 
cal for different cr'si^. Furthermore, it is apparent that 
the AH{M, AM) curves obtained from numerical simu- 
lations are asymmetric, in particular for small a. They 
show much larger deviations from the mean-field approx- 
imation on the left hand side, i.e. near negative satura- 
tion M = — 1. This can also be seen in the hysteresis 
loops themselves, where the curves seem to bundle up 
near negative saturation for small a values. 

In the right column of Fig. HI we show numerical val- 
ues for the deviation from redundancy measures (M) 
which are calculated from the simulated recoil curves 
shown in the left column. Due to the definition range of 
r,;j(M), only the recoil curve pairs — (1,2), (1,3), 
(1, 4) and (2, 3) produce data. From Fig.[3]^g), we see that 
for small a the deviation from data redundancy is quite 
substantial for the whole M definition range and for all 
the recoil curve pairs. For intermediate a, the deviation 
becomes smaller but is still visible, as shown in Fig.^h). 
For high a, the deviation is almost negligible in the whole 
M definition range and for all recoil curve pairs, as one 
can see from Fig. UKi) . Small non- vanishing values occur- 
ring at the boundaries of the M definition range are just 
due to the incompletely canceled finite size inaccuracies, 
as discussed in Sec. IIVBI 

For the other D{Hs) distributions, we observe very 
similar results as shown in Figs. EHZl Thus, one has to 
realize that for small cr, i.e. strong hysteron interactions, 
the AH{M, AM) method is not accurate. The mean- 
field approximation doesn't match the numerical result 
and deviations from redundancy are large. This is the ex- 
pected result because once coupling dominates the mag- 
netization reversal the mean-field approximation will not 
be valid any more. On the other hand, for large a, i.e. 
weak hysteron interactions, the AH{M, AM) method 
works very well, which is indicated by both the small 
deviation from redundancy and the match of the mean- 
field approximation to the numerical results. 



B. Emergent feature of the AH{M,AM) method 

The similarities in the failure of the AH{M,AM) 
method for different D{Hs) distributions indicates that 
this method may be not very sensitive to the particular 
type of distribution in general. To study this further, we 
plot the reliability measures against the tuning parame- 
ter a for all D{Hs) distributions in Fig. [51 

The fit quality measures Pd and are shown in 
Fig. [SKa) and (b), respectively. We see that Pd ap- 
proaches with increasing a, which means that for high 
cr, the input value of a can be recovered with very high ac- 
curacy by fitting the AH{M, AM) data to the mean-field 



fit function. In other words, the AH{M, AM) method 
works well for high a. Furthermore, we see that with 
increasing a, B? approaches 1 corroborating a successful 
fit of the AH{M, AM) data in this regime. The average 
deviation from redundancy (r) is shown in Fig. [Sfc). It 
is clearly seen that with increasing cr, r approaches 0, i.e. 
data redundancy is obtained, which is the key feature of 
the mean-field approximation. 



The reliability range of the AH{M, AM) method can 
be obtained from the reliability measures shown in Fig.[51 
We find that for all the four D{Hs) distributions, with 
the definitions of cr given in Sec. M the AH{M, AM) 
method works virtually perfect for cr > ctq with ctq being 
approximately equal to 20. Here, erg is just a rough cri- 
terion, above which the reliability measures have merged 
into their mean-field approximation values. It should be 
emphasized that the ctq values obtained from all three 
reliability measures are fairly consistent. If there are dif- 
ferences at all, r appears to show the highest sensitivity 
to deviations from the mean-field approximation, while 
seems to be slightly less sensitive. This is important, 
because r can be evaluated without any fit from a data 
set alone. So, the independent reliability test is the most 
sensitive measure and gives one confidence that not only 
the AH{M, AM) method fails in a well-defined way, but 
also that one can very reliably check for this failure mode. 



Finally, we note that the different r vs. a curves in 
Fig. ^c) for the different distribution types track each 
other almost exactly. However, this particular observa- 
tion dependents somewhat on how exactly we define a 
in the various D{Hs) distributions, because only for the 
Gaussian and Lognormal distribution are we using the 
natural definition given by the standard deviation. Thus, 
the curve collapse seen in Fig. [SKc) might be partially ar- 
tificial. Pd and R^ on the other hand do not exhibit such 
a collapse, not even for the Gaussian and Lognormal dis- 
tributions as is apparent from Fig. [5]^ a) and (b). The 
fact that the two kinds of reliability measures show dif- 
ferent behavior can be understood in the following way. 
As mentioned in the end of Sec. IIV Al and Sec. lIVBi fi- 
nite size inaccuracies at the definition range boundaries 
of the recoil curves will not affect the calculation of r 
very much due to the cancellation effect. But it will af- 
fect the calculation of the fit quality measures, both Pd 
and R^. Generally speaking, the shape of the hysteresis 
loops and recoil curves depends on the particular type 
of the chosen D{Hs). Consequently, the finite size in- 
accuracies will also depend on D{Hs)- As a result, we 
see different (similar) behaviors of the fit quality mea- 
sures (deviation from redundancy measure) for different 
D{Hs) distributions at small cr. In this sense, it is nat- 
ural to choose the deviation from redundancy measure r 
as the best measure to determine the reliability range of 
the AH{M, AM) method. 
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FIG. 4: Numerical results using a Gaussian distribution RSFHM. Rows: (Top, fig.a,d,g): a = 1.6. (Middle, fig.b.e.h); (7 = 5. 
(Bottom, fig.c,f,i): a = 50. Columns: (Left, fig.a,b,c): M{H) curves, main loop and 5 recoil curves. (Middle, fig.d,e,f): 
AH{M,AAI) curves for the 5 recoil curves: (solid lines) numerical result; (dotted lines) mean-field approximation. (Right, 
fig.g,h,i): M-dependent deviation from redundancy (rij(M)) for all the possible recoil curve pairs. 



VI. SUMMARY 

We study the AH{M,AM) method and its relia- 
bility by means of numerical simulations of the zero- 
temperature random switching field hysteron model. We 
present strong evidence that the AH{M, AM) method, 
which is based on the mean-field approximation, has a 
well-defined reliability range. This reliability range can 
be checked with two types of independent measures: de- 
viation from redundancy and fit quality. The former is 
the superior tool because it is calculated from the data 
set alone and is independent from any inaccuracies that 
might be induced by data fitting procedures itself. 
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APPENDIX A: MAPPING RSFHM TO RFIM 

For the zero-temperature non-equilibrium RFIM, a lo- 
cal metastable dynamics has been introduced by Sethna 
et alr^ to study the disorder-induced phase transition 
in the hysteretic behavior at T = 0: Initially, all spins 
point DOWN, as H is slowly increased from —oo to oo 
and decreased back to — (X), each spin flips deterministi- 
cally when its effective local field 

hf = jJ2sj+H + h, (Al) 
j 

changes sign. For the zero-temperature non-equilibrium 
RSFHM, we introduce a similar local metastable dynam- 
ics: Initially, all spins point DOWN, as H is slowly in- 
creased from —oo to oo and decreased back to — oo, each 
hysteron flips deterministically when its effective local 
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FIG. 5: Numerical results using a Lorentzian distribution RSFHM. Rows: (Top, fig.a,d,g' 
(Bottom, fig.c,f,i): a = 50. Columns: (Left, fig.a,b,c): M{H) curves, main loop and 5 
very large ho/<J ratio (~ 2 x 10^ for this system size 1000^, chosen to avoid negative Hs 
loop and all the recoil curves are extremely hard to see from the M{H) plot itself. The 
with log scale as shown in the middle column. (Middle, fig.d,e,f): AH{M,AM) curves 
numerical result; (dotted lines) mean-field approximation. Here, we see that normalized 
to the Gaussian D{Hs) case due to the large ho/a ratio. (Right, fig.g,h,i): M-dependent 
for all the possible recoil curve pairs. 



: a = 1.6. (Middle, fig.b.e.h): a = 5. 
recoil curves. Note that due to the 
), the differences between the major 
small differences will be more clear 
for the 5 recoil curves: (solid lines) 
AH values are very small compared 
deviation from redundancy {rij{M)) 



field 

Hf ^jY,Sj+H + sgn{Si) Hs, (A2) 

j 

changes sign. 

Considering the only slight difference between the 
RFIM and the RSFHM, it is possible to introduce a sim- 
ple mapping: hi <~* sgn(S'i) Hsi with Hsi > 0. This 
enables us to calculate the M{H) curve of the interact- 
ing hysterons by means of simple software adaptation. 



APPENDIX B: PROOF OF DATA 
REDUNDANCY IN MEAN-FIELD THEORY 



Ma = Mb ^ M (Bla) 
Mp = Mq = M - AMj + AM, (Bib) 
Afu = Mv = M - AMj (Blc) 

and M is an arbitrary value within [—1, 1 — AM^]. The 
hysteron distributions for the six states are shown in 
Fig. [51 Shadowed areas denote DOWN-hysterons {Si = 
— 1), while open areas indicate UP-hysterons {Si = -1-1). 
In the following, we use the symbol and for the 
density of DOWN-hysteron and UP-hysteron, respec- 
tively. For the states chosen from the upper branch of 
the major loop with magnetization Mq, it is easy to get 
n| = (1 — Mq)/2. For example, 



As shown in Fig. [31 we choose six states (A, B, P, Q, 
U and V) from the major loop and recoil curves with 



nx(^) = (l-M)/2 (B2a) 
n| (P) = (1 - M - AMi + AMj)/ 2 (B2b) 
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FIG. 6: Numerical results using a truncated Lorentzian distribution RSFHM. Rows: (Top, fig.a,d,g): a = 1.6. (Middle, 
fig.b.e.h): a = 5. (Bottom, fig.c,f,i): a — 50. Columns: (Left, fig.a,b,c): M(H) curves, main loop and 5 recoil curves. (Middle, 
fig.d,e,f): AH{M, AM) curves for the 5 recoil curves: (solid lines) numerical result; (dotted lines) mean-field approximation. 
(Right, fig.g,h,i): M-dependent deviation from redundancy {rij(M)) for all the possible recoil curve pairs. 



For states chosen from the recoil curves with magneti- It follows that 
zation Mq, there are two shadowed areas, which can be 
denoted as n^i and n\^2- Note that the left one n|i is just 
due to the distance from saturation at the starting point 
of the recoil curve, i.e. AM, so that = AM/2. For 
example, B and U are picked from the i-th recoil curve 
with distance from saturation AMi, so that 



"ii(S)="ii(C^) = AAf,/2 
Similarly, for Q and V, we have 

n^iiQ)^n^i{V)^AM,/2 



(B3) 



(B4) 



Combining this with the number conservation equation 
— '^ii + ^*i2 = (1 ~ Mq)/2, we find 



n^{A)=rn2{V) 
ni2{B) = ni2{Q) 



(B6a) 
(B6b) 
(B6c) 



ni2{B) 


= (1 


- M - AMi)/ 2 


(B5a) 


ni2{Q) 


-(1 


- M - AMi)/2 


(B5b) 


ni2iU) 


= (1 


- M - AM, + AMj)/2 


(B5c) 




= (1 


- M) /2 


(B5d) 



and 
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FIG. 7: Numerical results using a Lognormal distribution RSFHM. Rows: (Top, fig.a,d,g): a — 1.6. (Middle, fig.b.e.h): a — 5. 
(Bottom, fig.c,f,i): a = 50. Columns: (Left, fig.a,b,c): M{H) curves, main loop and 5 recoil curves. (Middle, fig.d,e,f): 
Aff(M, AA/) curves for the 5 recoil curves: (solid lines) numerical result; (dotted lines) mean-field approximation. (Right, 
fig.g,h,i): M-dependent deviation from redundancy (rij(M)) for all the possible recoil curve pairs. 



Fa + -ffi„t(M) = + ffi„t(M - AM,) 
iJp + i?,nt(M - A7\/j + AM,) =Hu + Hi^t{M - AM,) 

Hb + ffi„t(M) =Hq + iJint(M - AM, + AM,) 



(B7a) 
(B7b) 
(B7c) 



So, overall we find 



Q.E.D 



(Hb - Ha) + (Hy - Hu) - (iJq - Hp) = 
I 



(B8) 



APPENDIX C: A{H,AM) FORMULA FOR THE 
TRUNCATED LORENTZIAN DISTRIBUTION 



In numerical simulation, the random switching fields 
with any distributions are generated by a random number 
generator. To avoid any negative tails in the Lorentzian 
distribution, we can artificially suppress any negative 
random numbers and instead create another random 



number for the switching field until we get a positive 
one. The corresponding switching field distribution is 
then represented by the truncated Lorentzian distribu- 
tion Dl^{Hs): 



DuiHs) 





2 10 



C 



TV w'-^+MHs-hoY-' 

with C given by the 



for Hs < 

for iJs > '■^^^ 

condition that 
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FIG. 8: Reliability measures for the different D{Hs)'- (a) Pd- 
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C 



Specifically, one finds 



(C2) 



D(-Hs) ' 







FIG. 9: Hysteron distributions for the various states (A, B, P, 
Q, U and V) shown in Fig. |3l The shadowed areas symbolize 
DOWN-hysterons {Si = — 1), the open areas correspondingly 
UP-hysterons {Si = +1). 



For the truncated Lorentzian distribution, we define 
the disorder parameter a = w to make it comparable 
with the Lorentzian distribution. 

The integral function is given by: 



J — oo 

_i / 2{x - ho) 



C 

IT 



tan 



w 



so that 



From this result, we derive 



(C3) 



(C4) 



AHu{M,AM) = |{tan[(| + 0) 



1 - M 



-tan[(- + 0)( 



(C5) 



with = tan-i(^) 



as the analytic mean-field solution. We notice that due 
to the truncation induced asymmetry in the distribution 
itself, AHi,^{M, AM) depends on both w and Hq. 
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